Optical spectral variability of PG QSOs 
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Abstract. In previous studies we have shown that the 
optical variability of quasars increases, on average, with 
redshift. We explained this dependance in terms of a hard- 
ening of the spectrum during bright phases, coupled with 
the increase of the rest-frame frequency for increasing red- 
shift. We re-analize now these correlations on the basis of 
new light curves of PG quasars, recently published by the 
Wise Observatory group. 



1. Introduction 

Although variability plays a key role in constraining the 
models of the central engine of QSOs, little is known 
about the physical origin of luminosity variations. The 
most diverse variability mechanisms have been proposed 
in the past, including supernovae explosions (Aretxaga et 
al. 1997), instabilities in the accretion disk (Kawaguchi 
et al. 1998), and gravitational lensing due to intervening 
matter (Hawkins 1993). So far, most of the information 
about the characteristics of variability derives from the 
statistical analysis of single band light curves of magni- 
tude limited QSO samples. The correlation of variability 
with either intrinsic luminosity (v — L) or redshift (v — z) 
is affected by the strong correlation between luminosity 
and redshift (L — z), present in these samples. The results 
of these analyses also depend on the specific variability 
index adopted, as shown by Giallongo et al. (1991), who 
found a positive v-z correlation through a variability index 
based on the rest-frame structure function, later confirmed 
by Cristiani et al. (1996). The increase of variability with 
the rest-frame observing frequency found by Di Clemente 
et al. (1996), supports the suggestion by Giallongo et al. 
(1991) that QSOs at high redshift appear more variable, 
on average, since they are observed in a higher rest frame 
frequency. The dependence of variability on frequency is 
associated with the hardening of the spectral energy distri- 
bution (SED) during the bright phases (Cutri et al. 1985, 
Kinney et al. 1991, Edelson et al. 1990). Trevese et al. 
(1999) have shown that a hardening of the SED in the 
bright phase occurs, on average, in the statistical, magni- 
tude limited, sample of AGNs of the SA 57 (Trevese et al. 



1989, 1994, Bershady et al. 1998). They also found that 
the slope a and its variations associated with the luminos- 
ity changes, are consistent on average with temperature 
changes of a black body. 

1.1. Data and analysis 

We present preliminary results of a new analysis of the 
spectral variability of QSOs in the optical band. It is based 
on a data set made available by Giveon et al. 1999, con- 
sisting of the light curves of a sample of 42 nearby and 
bright QSOs (z < 0.4, B < 16 mag) belonging to the 
Palomar-Green (PG) sample, which were monitored for 7 
years, with a typical sampling interval of 39 days, in the 
Johnson Cousins B and R bands, with the 1 m telescope of 
the Wise Observatory. Here the quantities adopted for the 
analysis are defined. The results are discussed in Section 
2, together with the relevant figures. 

1.2. Variability index 

We define the structure function as (see Di Clemente et 
al. 1996) Si = [7r/2((|m(i + r) - m(t)\) 2 - cr 2 )] 1 / 2 , where 
m(t) is either the B or the R magnitude, t is the rest- 
frame time, t is the time lag between the observations, 
o~ n is the relevant r.m.s. noise and the angular brackets 
indicate the average taken over all the pairs of observa- 
tions separated by a time interval r ± At. The value of 
At is the result of a trade-off between time resolution and 
statistical uncertainty. In the following we define, for each 
QSO, two variability indexes vb and Vr equal to the 
relevant structure function in the B and R bands, com- 
puted for r = 2 ± 0.6 yr. The average values of vr and vb 
over the ensemble of 42 QSOs, versus the relevant average 
rest-frame frequency, are reported in Figure 1, which is 
adapted from Di Clemente et al. (1996). 

1.3. Instantaneous spectral slope 

For each QSO we compute the instantaneous spectral 

slope a(t) = log(/ 1/B // I/J j)/log(i/ B /i'H) ) regarding as syn- 
chronous the observations within a time lag of 9 hours. 
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Emission line corrections are not considered in the present 

. , . 0.4 

preliminary analysis. 

1.4- Spectral variability parameter 

Considering all possible pairs of observations, separated by 3 
a time interval = Tj — Tj , we define a spectral variability 
indicator: (3{t) = [a(t+r) -a(t)}/ [log L B (t+r) -log L B (t)} 
where Ls(t) is the absolute luminosity in the blue band. 
In the present analysis we consider, for each object, the 
spectral variability parameter (3 m defined as the av- 
erage of all the (3(t) values for r < 1000 days. 



2. Results 

2.1. Variability- frequency correlation (v — v) 

Fig. 1 shows the variability indexes corresponding to var- 
ious frequencies. The stars represent the indexes Vr and 
vb computed in the present analysis. The other points, 
taken from Di Clemente et al. 1996, derive from different 
statistical QSO samples as indicated. The new points are 
consistent with the general trend and confirm the value of 
the slope dSi/dlogv reat = dS 1 /dlog(l + z) = 0.25-0.30 
which accounts for the v — z correlation found by Gial- 
longo et al. (1991). A v — z correlation is not detectable 
for the whole sample (Giveon et al. 1999), due to the 
small redshift range (z < 0.4). If we restrict to an ab- 
solute magnitude bin —23.5 < Mb < —22.5, to reduce 
the combined effect of the (positive) L — z and (nega- 
tive) v — L correlations, which attenuates the apparent 
the v — z correlation, we find a partial correlation coeffi- 
cient r V)Z = 0.39, marginally significant [P(> r) = 0.09], 
despite the small number (19) of QSOs in the selected 
magnitude bin. Therefore, we can conclude that the in- 
trinsic v — z correlation is present, and the trend shown in 
Fig. 1 confirms its interpretation in terms of an increase 
of variability for increasing rest-frame frequency. 



2.2. Spectral slope - luminosity (a — Lb) correlation 

In Fig. 2 the instantaneous slope a(t) of each QSO is plot- 
ted as a function of the relevant blue luminosity L b ■ The 
regression lines a vs. Lb are also reported for each QSO. 
Two features appear in the figure: (i) Intra-QSO corre- 
lation (Aa — A log L) , i.e. a trend (with a few exceptions) 
of increasing a for increasing log Lb appears, indicating 
that QSO SEDs are harder in the bright phase (and vice 
versa); (ii) Inter-QSO correlation (a — logL), i.e. also 
the average a of each QSO is positively correlated with 
its average luminosity. 

These two effects suggest that both the change of a 
during luminosity variations and the average increase of a 
from faint to bright QSOs have the same physical origin: 
the increase of the temperature of the emitting gas. As 
an example the a(t) versus log Lb (i) lines computed for 
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Fig. 1. Variability indexes as a function 
quency adapted from Di Clemente et al. 
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Fig. 2. Instantaneous slope a(t) of each QSO, plotted as 
a function of the relevant blue luminosity L b ■ The regres- 
sion lines a vs. logL^ are also reported for each QSO. 



three emitting black bodies of different areas are reported 
in Fig. 2. The temperature increases along each line from 
bottom left to top right. 
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Fig. 3. Spectral variability parameter versus spectral 
slope. 

2. 3. Spectral variability versus spectral slope (3 — a 

Fig. 3 shows the spectral variability parameter (3 m of each 
QSO, versus its average spectral slope a m . The continuous 
curve represents f3 versus a for a family of black bodies of 
different temperatures, increasing from left to right. This 
is not a fit to the data points, since there are no free pa- 
rameters in the curve. Comparison with the data suggests 
that the spectral slope and the slope variations are roughly 
consistent with a black body of tempeature T = 2.5 • 10 4 
K, subject to small temperature changes responsible of 
both the luminosity and spectral slope changes. 



— the two latter points suggest that the gas temperature 
is responsible for both the effects; 

— temperature variations of a black body with an average 
temperature of the order of 2.5 • 10 4 K are consistent 
for the intra-QSO a(t) — Lb if) correlation. 

— different average black body temperatures may ac- 
count for the inter-QSO a — Lb correlation. 
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3. Conclusions 

The analysis of the spectral variability of the present sta- 
tistical subsample of PG QSOs shows that: 

— the amplitude of variability increases with the rest- 
frame frequency; 

— a (positive) v — z correlation appears in the present 
sample, despite the (negative) v — L and the (positive) 
L — z correlations; 

— the increase of variability with frequency is consistent 
with the trend previously found by Di Clemente et al. 
(1996) on the basis of different QSO samples and con- 
firms the explanation of the v — z correlation suggested 
by Giallongo et al. (1991); 

— the spectral slope increases on average in the bright 
phases (intra-QSO a — L correlation); 

— the spectral slope is steeper, on average, for fainter 
QSOs (inter-QSO a — L correlation); 



